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The Sine quadrature and the Sine indefinite integration are approximation formulas for 
definite integration and indefinite integration, respectively, which can be applied on any 
interval by using an appropriate variable transformation. Their convergence rates have been 
analyzed on any interval. In addition, with the aim of verified computation, more explicit 
| error bounds that are computable have been recently given on the finite interval. In this 

paper, such explicit error bounds are given in the remaining case: on the semi-infinite interval 
and the infinite interval. 
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1 Introduction 



, The "Sine quadrature" is an approximation formula for the integral over the whole real line, 

expressed as 

,oo N 

/ F(x)dx^h J2 F(kh), (1.1) 
J -°° k=-M 

where M, N, h are selected appropriately depending on n. This approximation is also called 
the (truncated) "trapezoidal formula." It is well known that the formula (|1.1|) can achieve quite 
a fast, exponential convergence. Furthermore, its optimality is proved under a certain class of 
functions [4,6j. Here, there are two important conditions to be satisfied: (i) the interval of 
^ integration is (— oo, oo), and (ii) |-F(x)| decays exponentially as x — > ±oo. In other cases, users 

should employ an appropriate variable transformation t = ip(x), i.e., the given integral should 
be transformed as f(t) dt = Jf^ f(ip(x))'i^'(x) dx, so that those two conditions are met. 
In this regard, Stenger [S] considered the following four typical cases: 

1. (a, b) = (— oo, oo), and \f(x)\ decays algebraically as x — > ±oo, 

2. (a, b) = (0, oo), and \f(x)\ decays algebraically as x — > oo, 

3. (a, b) = (0, oo), and decays (already) exponentially as x — > oo, 

4. The interval (a, b) is finite, 
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and gave the concrete transformations to be employed in all cases: 



ip SE i(t) = sinht, 

ip SE2 (t) = e*, 

ipSE 3 (t) = arcsinh(e*), 

b — a ( t\ b + a 
VsE4(t) = — — tanh I — I H — , 



which are called the "Single-Exponential (SE) transformations." Takahasi-Mori [TJ proposed 
the following improved transformations: 

ip um (t) = sinh[(7t/2)sinh£], 
VW*)=e (7t/2)sinhi , 

, . . b — a f 71 . \ b + a 
VdeS) = — tanh I — sinhi 1 H — , 

which are called the "Double-Exponential (DE) transformations." In addition, in the case 3, 
another DE transformation 

^ DE3 (t)=log(l+e( n / 2 ) sinh4 ) 

was proposed [2] so that its inverse function can be explicitly written with elementary functions 
(whereas iposs^t) cannot). 

Error analyses of the Sine quadrature combined with ip S Ei(t), • • • , ips™ (t) have been given [S] 
in the following form: 

|Error| < Ce"^, 
and for i/> DE i(t), ipnmit), ipoE3(t), tpom(t), their error analyses have been given [S] as 

|Error| < C &~ 2ndnl lo s( M ™/^) , (1.2) 

and for V ; DE3t(^)) a ls° given [5] as 

|Error| < C e~ 2ndn / W^/n) , (1.3) 

where [i indicates the decay rate of the integrand, and d denotes the width of the domain in 
which the transformed integrand is analytic, and C is a constant independent of n. In view of 
the inequalities above, we notice that the accuracy of the approximation can be guaranteed if 
the constant C is explicitly given in a computable form. In fact, the explicit form of C has been 
revealed in the case 4 (the interval is finite) [3], and the result was used for verified automatic 
integration [S]. 

The main objective of this study is to reveal the explicit form of C's in the remaining cases: 
1-3 (the interval is not finite), which enables us to bound the errors by computable terms. 
As a second objective, this paper improves the DE transformation in the case 3. Instead of 

^ DB3 (i) or ^DB 3 f(i), 

^ DE3t (t) = log(l+e" sinht ) 
is introduced in this paper, and it is shown that the error is estimated as 

|Error| < C e- 2ndn ' Xo ^ Adn ^\ (1.4) 
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with clarifying the constant C. The rate of (|1.4p is better than (|1.2p and f j 1 . 3 j> . Furthermore, by 
using Sugihara's nonexistence theorem [BJ, it can be shown that ^ DE3 } is the best among possible 
variable transformations in the case 3 (although the point is not discussed in this paper and left 
for another occasion). 

In addition to the "Sine quadrature" described above, similar results can be given for the 
"Sine indefinite integration" for indefinite integrals f(t) dt, which is also examined in this 
study. 

The remainder of this paper is organized as follows. Main results of this paper are stated in 
Sections [2] and [3j new error estimates for the Sine quadrature are presented in Section [21 and 
for the Sine indefinite integration in Section [31 All proofs of the presented theorems are given 
in Section [5j Numerical examples are shown in Section |H 



2 Error Estimates with Explicit Constants for the Sine Quadra- 
ture 

In this section, new error estimates for the Sine quadrature are stated comparing with the 
existing results. Let us introduce some notation here. Let S>d be a strip domain defined by 
S! d = {C G C : |ImC| < d} for d > 0. Furthermore, let f d " = ^n{(eC: ReC < 0}, 
= &d n {C G C : Re ( > 0}. In all theorems in Sections [2] and [3J d is supposed to be a 
positive constant with d < n/2. For a variable transformation ip, ip{@d) denotes the image of 
3> d by the map ip, i.e., ip(@d) = {z = -0(C) : C G &d}- Let I\ = (— oo, oo), I 2 = h = (0, oo), and 
let us define the following three functions: 

E 2 (z;a,(3) 

We write Ei(z;j,j) as 7) for short. 



1 

(l + Z 2)( 7 +l)/2' 

z a-l 
(1 +z 2^(a+P)/2- 
a-l 

l + z 



2.1 Existing and New Error Estimates for the Sine Quadrature with the SE 
Transformation 

Existing error analyses for the Sine quadrature with V'seij ^se2, V'sbs are written in the following 
form (Theorems 12.11 and I2.2[) . 

Theorem 2.1 (Stenger Theorem 4.2.6]). Assume that / is analytic in ipsEi(^d), and there 
exist positive constants K, a, and j3 such that 

|/(*)| <K\E x {z;a)\ (2.1) 



for all z G ip SE1 (£? d ), and 

\f(z)\ <K\E^z;P)\ (2.2) 
for all z G tp SE i(^). Let fi = min{a, /?}, let h be defined as 

1 2nd , , 

h = J , (2.3) 
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and let M and N be defined as 



'M = n, N=\an/(3] (if \i = a), 
= n, M = \(3n/a] (if fx = (3). 



(2.4) 



Then there exists a constant Ci, independent of n, such that 



N 



h 



f(t)dt~h J2 fU>BBl(kh)W am (kh) 



k=-M 



(2.5) 



Theorem 2.2 (Stenger [SJ Theorem 4.2.6]). The following is true for i = 2, 3. Assume that / 
is analytic in tpsEi(^d), an d there exist positive constants K, a, and /3 such that 



|/(z)| <K|^(z;a,/3)| 



(2.6) 



for all z G V'sEi(-^rf)- Let /x = min{a, /3}, let /i be defined as (|2.3|) . and let M and N be defined 
as ()2.4p . Then there exists a constant Ci, independent of n, such that 



N 



f(t)dt-h J] f& SEi (khW SEi (kh) 



k=-M 



(2.7) 



This paper explicitly estimates the constant Cj's in (|2.5|) and (|2.7|) as follows. 



Theorem 2.3. Let the assumptions in Theorem EH] be fulfilled. Furthermore, let v = max{a, /?}. 
Then the inequality (|2.5|) holds with 



Ci 



2 V+1 K 



fl 1 (1 - e -V / 2^){cOs4' 



+ 1 • 



Theorem 2.4. Let the assumptions in Theorem l2.2l be fulfilled. Then the inequality (|2.7|) holds 
with 



A* )(l _ e - v / 2^){cos ci}^)/ 2 



Ai [ (i _ e -V^){ CO sdY a +P)/ 2 
where c a 4 is a constant defined by 



+ 1 



+ 1 



! \(l-a)/2 



1 



(if < a < 1) 
(if 1 < a). 
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2.2 Existing and New Error Estimates for the Sine Quadrature with the DE 
Transformation 

Existing error analyses for the Sine quadrature with V'dei, V ; de2, i^DE3 and ^ D E3t are written in 
the following form (Theorems 12.51 and 12. 6|) . 
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Theorem 2.5 (Tanaka et al. [81 Theorem 3.1]). The following is true for i = 1,2, 3. Assume 
that / is analytic in ^DEi(^d); an d there exist positive constants K and fj, (with \i < 1 in the 
case i = 3) such that 

\f(z)\<K\Ei(z;ii)\ 

for all z £ ipuE t (^d)- Then there exists a constant C, independent of n, such that 



where 



f(t)dt-h J2 f(^(kh))i;' DE1 (kh) 

k=—n 

log(8cin//j) 



^ (j ^—2ndn/ log(8dn//j,) 



n 



(2.9) 



Theorem 2.6 (Tanaka et al. [SJ Theorem 3.1]). Assume that / is analytic in V'DE3t(^d)) an d 
there exist positive constants K and /U with /x < 1 such that 

|/(z)| < if |£ 3 (^)l 

for all z £ ^DE3t(^rf)- Then there exists a constant C, independent of n, such that 



J~3 



f(t)dt-h J2 f(^{kh)W^{kh) 

k=—n 



g— 2ndn/ log(2ndn / /i) 



where /i = log(27tdn/yu)/n. 

Remark 2.1. As for Theorem 12.51 with i = 3 and Theorem 12.61 although the condition < 1' 
is not assumed (only \i > is assumed) in the original paper [8J, that condition is necessary to 
avoid the case where \z/(l + z)^' 1 = oo at z = — 1 (see E%(z; fx)). 

As for the case 1 (Theorem 12.51 with i = 1) and the case 2 (Theorem 12.51 with i = 2), this 
paper not only explicitly estimates the constant C's, but also generalizes the approximation 



formula from J2k=-n *° S 



A" 



-Af 



as stated below. Here, x 7 is defined for 7 > by 



arcsinh 



2~TVy 



arcsinh(l) 



(if < 7 < l/(27t)), 
(if 1/(2tt) < 7). 



Theorem 2.7. Assume that / is analytic in ipv>Ei{^d)i an d there exist positive constants K, 
a, and f3 such that (|2.ip holds for all z G ^dei(^), and (|2.2|) holds for all z € ip OE1 (@£). Let 
/i = min{a, /?}, let f = max{a, /3}, let /i be defined as (|2.9p . and let M and A be defined as 



(2.10) 



'M = n, N = n- [log(P/a)/h\ (if/z = a), 
AT = n, M = n- \log{a/P)/h\ (if /J = /?)• 



Furthermore, let n be taken sufficiently large so that n > (i/e)/(8d), Mh > x a , and iV/i > xp 
hold. Then it holds that 



N 



f(t)dt-h £ /(VW^M E1 (^) 



where Ci is a constant independent of n, expressed as 



< e -27tdn/ log(8dn/V) 



2 y+1 if 



(1- e-^ e / 4 ){cos(f sin d)}" cos d 



+ e 
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Theorem 2.8. Assume that / is analytic in '0 D E2(^d) I and there exist positive constants K, 
a, and (3 such that (|2.6|) holds with i = 2 for all z E "0DE2(^d)- Let ti = min{a, let 
v = max{a, /?}, let h be defined as (|2.9|) . and let M and N be defined as (|2,1U|) , Furthermore, 
let n be taken sufficiently large so that n > (i>e)/(8d), Mh > x a , and Nh > X/j hold. Then it 
holds that 

N 

f(t)dt-h J2 f(^(kh))^ 2 (kh) 

h k=-M 

where C2 is a constant independent of n, expressed as 



< Q e -2ndn/ log(8dn/fi) ^ 



c = 2 Ji J 2 /4 1 

2 /i \ (1 - e-^ e / 4 ){cos(f sin d)}(<*+P)/ 2 cos d J ' 

As for the case 3 (Theorem 12 . 51 with i = 3 and Theorem l2.6p . this paper employs the improved 
variable transformation ip OE3i as described in the introduction. Then this paper gives the error 
estimates in the similar form to Theorems 12.71 and 12.81 

Theorem 2.9. Assume that / is analytic in V J DE3t(^d)j an d there exist positive constants K, 
[3 and a with a < 1 such that (|2.6|) holds with i = 3 for all 2 G VtoEstO^d)- Let ti = min{a, /?}, 
let = max{a, /?}, let /i be defined as 

fc= Iog(4*i/ M ) (2 . n) 
n 

and let M and N be defined as ()2.10p . Furthermore, let n be taken sufficiently large so that 
n > (ye) /{Ad), Mh > x a , and Nh > hold. Then it holds that 



N 

f(t)dt-h ]T f(i; DE3t (kh))^, t (kh) 
13 k=-M 

where Cq is a constant independent of n, expressed as 



< (J e -27rdn/ log(4dn/^t) 



~ M-a 



£. = ^ I 2 (c~d) 7t(l-a+6«/)/12 

3t /i \(1- e-^ e / 2 ){cos(f sin d)} a +0 cos d j 
and where is a constant expressed by using c<i = 1 + {l/cos(^ sind)} as 

. 1 + lQg(l + Cd) 

Cd = — ; — 77— 7— c rf . (2.12) 

log(l + Cd) 

3 Error Estimates with Explicit Constants for the Sine Indefi- 
nite Integration 



Sine indefinite integration is an approximation formula for the indefinite integral [T], expressed 



as 

N 



k=-M 

Here, the basis function J(k, h) is defined by 



j F(x)dx^ F(kh)J(k,h)(£), ^1. (3.1) 



J(k, h)(x) = h\\ + - Si[n(x/h - k)] 

{ 2 7T 



where Si(x) is the so-called sine integral, defined by Si(x) = {sin(a) / a} da. The approxima- 
tion (|3.ip can be combined with the SE transformation or the DE transformation [2] similar 
to the Sine quadrature (|1.1|) , This section is devoted to presenting the error estimates for the 
formulas. 
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3.1 New Error Estimates for the Sine Indefinite Integration with the SE 
Transformation 

This paper gives new error estimates for the Sine indefinite integration with ip SE1 , ip S E2, and ip SE3 
in the following form (Theorems 13.11 and 13. 2[) . 

Theorem 3.1. Assume that / is analytic in ipsm{@d)i arid there exist positive constants K, 
a, and (3 such that ([2TT]) holds for all z G ip 3m (@~[), and (f!T2]) holds for all z G ifj SE1 (^). Let 
\x = min{a, /?}, let v = maxja, /?}, let h be defined as 



y /in' 

and let M and N be defined as ([53]). Then it holds that 



(3.2) 



sup 



N 



f(t)dt- f(A E i(kh))^ E1 (kh)J(k,h)(^- E \(r)) 



' k=-M 

where C\ is a constant independent of n, expressed as 



< Cie 



- \J nd/in 



Ci 



2 V+1 K 



71 



^ _ e -2V^M){cos4" V d/i 



+ 1.1 



Theorem 3.2. The following is true for i = 2, 3. Assume that / is analytic in V'sEi(^d)) and 
there exist positive constants K, a, and j3 such that (|2.6p holds for all z G ipsEi(^d)- Let 
/i = min{a, /?}, let /i be defined as (|3.2|) . and let M and iV be defined as (|2.4|) . Then it holds 
that 

AT 

/(t)dt- 2 /(^s El (A:/ i ))^E l (^)^»(^"E 1 I W) 
k=-M 



sup 

TGl; 



-■^/nd/in 



where C2 and C3 are constants independent of n, expressed as 



Co 



2K 



n 



H e - 2 V^){cos d}( a +P)/ 2 V d ^ 



2K 



2 (a+ y 8-l)/2 



7T 



M [ (1 _ e _2 v^){cos d}( a +/ ? )/ 2 V d/i 
and where c a d is a constant defined in (12.81) . 



+ 1-1 > , 



+ 1-1 > , 



3.2 New Error Estimates for the Sine Indefinite Integration with the DE 
Transformation 

This paper gives new error estimates for the Sine indefinite integration with V'dei) V'de2> and 
^DEat in the following form (Theorems E3H32} • Let = [e~ ndn / lo ^ dn / ^ \og{Adn/ n)]/n for 
short. 

Theorem 3.3. Assume that / is analytic in V'dei(^)) and there exist positive constants K, 
a, and j3 such that (|2.1|) holds for all z G ^ DE i(^"), and (|2.2p holds for all z G ^ DE1 (^j"). Let 
\x = minja, /?}, let 1/ = maxja, /3}, let /i be defined as (|2.11|) . and let M and N be defined 
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as ()2.10p . Furthermore, let n be taken sufficiently large so that n > (ue)/(4d), Mh > x a , and 
Nh > X/3 hold. Then it holds that 



sup 

T€/l 



A? 



f(t)dt- f(^i(kh))iP' DE1 (kh)J(k,h)^-i(T)) 

k=-M 



< CieT 



where C\ is a constant independent of n, expressed as 
2 U+1 K f 1 



Ci 



fid 1(1- e- 7t ' ie / 2 ){cos(2 sin d)}^ cos d 



+ e 



7r(a+/3)/4 



Theorem 3.4. Assume that / is analytic in V'de2(^)) and there exist positive constants K, 
a, and (3 such that (|2.6|) holds with i = 2 for all z G i^wi^d)- Let ,u = min{a, let 
f = max{a, let /i be defined as (|2.11|) . and let M and A be defined as (|2.1U|) . Furthermore, 
let n be taken sufficiently large so that n > (i/e)/(4<i), M/i > x a , and A/i > xp hold. Then it 
holds that 



sup 

T&I 2 



N 



/(t)dt- £ f(^(kh))^(kh)J(k,h)(^(r)) 



k=-M 



where C2 is a constant independent of n, expressed as 
r 2A f l 



/xd 1 (1 - e"^ e / 2 ){cos(f sin d)}( a +^)/ 2 cos d 



+ e 



7t(a+/3)/4 



Theorem 3.5. Assume that / is analytic in ipnE3t(^d) , an d there exist positive constants K, f3, 
and a with a < 1 such that (|2.6p holds with i = 3 for all 2; £ V'des^^)- Let fi = min{a, (3}, let 
v = max{a, (3}, let h be defined as h = \og(2dn/ fi) /n, and let M and A be defined as ()2.10l) . 
Furthermore, let n be taken sufficiently large so that n > (ue)/(2d), Mh > x a , and Nh > xp 
hold. Then it holds that 



sup 

rG/ 3 



N 



f(t)dt- /(^DE3 t (A:/ l )X E 3 t (fe/ 1 )J(fc,/ l )(^ D : E 1 3 t (r)) 



-M 



where is a constant independent of n, expressed as 



r - 2K 



fid \ e-^ e ){cos(f sind)} a +^ cosd 
and where ca is a constant defined in (|2.12|) . 



+ e 



7r(l+5a+6/3)/12 



4 Numerical Examples 

In order to numerically confirm the results presented in Sections [5] and El let us consider the 
following three examples. 

Example 4.1 (Case 1 [2J). Consider the function f x (t) = V?,/{2n(t 2 + t + 1)) and its defi- 
nite/indefinite integral on I\: 



/oo 
/i(t)dt=l, 
-00 

£ o /i(*)d* = i + i«c tan {^(r + i)}. 



(4.1) 
(4.2) 
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The integrand f\ satisfies the assumptions in Theorem 12.31 with a = = 1, d = 3/4, and 
K = y/3e. In addition, fi satisfies the assumptions in Theorem 12.71 with a = /3 = 1, d = 7t/7, 
and K = 8V3/e. 

Example 4.2 (Case 2 [2]). Consider the function f2(t) = 2/(7t(l+t 2 )) and its definite/indefinite 
integral on I2: 

rod 

/ / 2 (t)dt = l, (4.3) 
Jo 

f T 2 

/ f 2 (t)dt = -arctan(r). (4.4) 
Jo 7t 

The integrand f% satisfies the assumptions in Theorem 12.41 (i = 2) with a = /3 = 1, d = cosh(l), 
and K = 2/n. In addition, f% satisfies the assumptions in Theorem 12 . 81 with a = f3 = 1, d = 3/2, 
and K = 2/n. 

Example 4.3 (Case 3 |7J). Consider the function fs(t) = e~^ 1+ ^ /(1+t) and its definite/indefinite 
integral on ^3: 

/ / 3 (i)dt = Ei(l), (4.5) 
Jo 

f T / 3 (t)dt = Ei(l)-r(0,l + r), (4.6) 

where Ei(a;) is the exponential integral, and T(s,x) is the incomplete gamma function. The 
integrand f$ satisfies the assumptions in Theorem 12.41 (i = 3) with a = /3 = 1, d = 3/2, and 
K = e~ l . In addition, / 3 satisfies the assumptions in Theorem 12.91 with a = /3 = 1, d = log(7t), 
and K = e. 

The numerical results are shown in Figures [TH3 All programs were written in C++ with 
double-precision floating-point arithmetic, and the GNU Scientific Library was used for comput- 
ing special functions. In Figure^ 'maximum error' denotes the maximum value among the errors 
investigated on the following 403 points: r = 0, +2" 100 , ±2"", . . . , +2" 1 , ±2°, +2 1 , . . . , ±2 100 . 
Similarly, in Figures 0] and El the errors were investigated on the 201 points (just the positive 
points of above), and its maximum is shown on the graphs. In each graph, we can see that the 
error estimate by the presented theorem (dotted line) surely bounds the actual error (solid line). 

5 Proofs 

5.1 In the Case of the SE Transformation 

Let us have a look at the sketch of the proof first, taking the Sine quadrature as an example. 
Let F(x) = ip SEi (x)ip' SEi (x) (recall that we employ the SE transformation t = ip SEi (x)). Then we 
have to evaluate the following term: 

N N 

/ f(t)dt-h ftt> SEi (khM Ei (kh) = / F{x)dx-h ]T F(kh) 

J Ii 1 »r J— OO 7 »r 



k=-M 



k=-M 



For the estimation, the function space defined below plays an important role. 

Definition 5.1. Let L, R, a, f3 be positive constants, and d is a constant with < d < n/2. 
Then L| E ^ Qj g(^) denotes a family of functions F that are analytic on f^, and satisfy for all 
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Figure 1. Error of the Sine quadrature Figure 2. Error of the Sine indefinite integra- 
tor (|4.ip and its estimate. tion for (|4.2|) and its estimate. 
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Figure 3. Error of the Sine quadrature Figure 4. Error of the Sine indefinite integra- 
tor (|4.3p and its estimate. tion for (|4.4|) and its estimate. 
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Figure 5. Error of the Sine quadrature Figure 6. Error of the Sine indefinite integra- 
tor (|4.5p and its estimate. tion for (|4.6|) and its estimate. 



10 



C e &d and i£l that 



\H0\ < 

\F(x)\ < 



L 



|1 + e -2C|a/2| 1 + e 2C|/3/2' 

R 

(1+ e- 2 *)<*/ 2 (l+ e 2x )?/ 2 ' 



(5.1) 
(5.2) 



If F belongs to this function space, the error of the Sine quadrature is estimated as below. 
The proof is omitted since it is quite similar to the existing theorem for the case 4 j3[ Theo- 
rem 2.6]. 



Theorem 5.1. Let F G L| E fi a p(S>d), let A 4 = min{a, /?}, let h be defined as 
and N be defined as (|2"H). Then it holds that 



, and let M 



/ F(x)dx-h F (kh) 
J -°° k=-M 



< 



2L 



(1- e-V^){ cos d}( a +P)/ 2 



+ R 



- \J 2nd/m 



This theorem states the desired error estimates with explicit constants for the Sine quadrature 
(if F 6 L| E r q si^d), which is not confirmed yet). 

For the Sine indefinite integration, the next theorem holds. The proofs is also omitted since 
it is quite similar to the existing theorem for the case 4 [31 Theorem 2.9]. 



Theorem 5.2. Let F £ L| E i?Q/3 (^), let /i = min{a, let h be defined as 
and N be defined as (12.41). Then it holds that 



sup 



and let M 



rt N 

/ F(x)dx- J2 F(kh)J(k,h)(C) 



< 



k=-M 

L 



(1 



-2y/ndn 



){cosd}( Q +/ 3 )/ 2 



+ l.LR 



^— nd/m 



This theorem states the desired error estimates with explicit constants for the Sine indefinite 
integration. 

In view of Theorems 15.11 and 15.21 our project is done by checking F £ ^LRafsi^d) in each 
case: 1, 2, 3. Let us check them one by one in what follows. 

5.1.1 Proofs in the Case 1 (Theorems [273] and [3TT]) 

The claims of Theorems 12.31 and 13.11 follow from the next lemma. 

Lemma 5.1. Let the assumptions in Theorem 12 .31 or Theorem [XT] be fulfilled. Then the function 
F(0 = /(V> 8E i(CMei(C) belongs to L s L \ a ^ d ) with L = 2 v K/{cos d}^~^ 2 and R = VK. 



Proof. Consider the case a < f3 first. From the inequality (|2.1|) . it follows that 

K 



|F(C)| <tf|£i(fei(C);a)lK E i(C)l 



|1 + e -2C|Q/2| 1 + e 2C|/3/2 

for £ € @a with Re£ < 0, and from the inequality (|2.2|) . it follows that 

K 



2 Q 1 + e 



2C|C8-q)/2 



|F(C)| <K|i?i(^ E1 (C);/?)ll^ B1 (C)l 



2^ 



|1 + e -2C|a/2|l + e 2C|/3/2 (]_ + e -2C|09-a)/2 
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for ( E &d with Re£ > 0. Setting £ = x + iy with x < 0, we have 
2*| l + e 2 ^"")/ 2 = 2 a (l + e 2 *)^-^/ 2 {l 



. 2 1 (Z 3 -")/ 4 

sin y 



cosh 2 x 

< 2 a (l+ e °)^- Q )/ 2 {l-0} (/3 - a)/4 

2^ 



= 2 (°+/3)/2 < 2 /3 < 

Furthermore, setting £ = x + iy with x > 0, we have 
2^ 2^ 



{cosy}^"")/ 2 ' 



|1 + e -2f|08-a)/2 



< 



(1 + e -2x)C9-a)/2 | 
2^ 



t sin* W ] (Z 3 -")/ 4 
cosh 2 z J 



2^ 



d + O^/aji-^L} 

Thus, since \i = a and v = /3 in this case, it holds for all £ £ £^ that 

ft" 2^ 

1^(01 < 



)/ 4 {cosy}^-")/ 2 



1+ e- 2 q°/ 2 |l+ e 2 ^/ 2 {cosd}^)/ 2 ' 



and it holds for all x S R that 
|F(x)| < 



K 



2 u 



(1 + e" 2:E ) a / 2 (l + e 2 *)^ 2 {cos 0}(^-^)/ 2 ' 
In the case a > (3, the same inequalities hold. This completes the proof. 



5.1.2 Proofs in the Case 2 (Theorems [273] and [372] with i = 2) 

The claims of Theorems 12.41 and 13.21 (i = 2) follow from the next lemma. 

Lemma 5.2. Let the assumptions in Theorem 12.41 or Theorem 13.21 be fulfilled with i = 2. Then 
the function F(C) = /(^se2(C))^se 2 (C) belongs to ^R^pi^d) with L = K and R = K. 

Proof. From the inequality (|2.6p with i = 2, immediately (|5.1|) and (|5.2p hold with L = R = 
K. U 



5.1.3 Proofs in the Case 3 (Theorems [23] and [3T2] with i = 3) 

The claims of Theorems 12.41 and 13.21 (i = 3) follow from the next lemma. 

Lemma 5.3. Let the assumptions in Theorem 12.41 or Theorem 13.21 be fulfilled with i = 3. Then 
the function F(Q = /W>seb(CMes(0 belongs to L| E flQ /3 (^) with L = 2^- v )l 2 c a4 K and 
R = K, where c a a is a constant defined in 



For the proof, let us prepare some useful inequalities (Lemmas I5.4H5.7|) . 



Lemma 5.4. For all £ G ^t/2) it holds that 



1 

72 



< 



arcsinh( e^ 



1 + e< 

Furthermore, for all x 6 R, it holds that 

arcsinh( e 



1 + arcsinh( e^) 



< y/2 


4 


1 + ef 





< 



1 + arcsinh( e x ) 1 + e x 



(5.3) 



(5.4) 
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Proof. Consider (|5.4|) first, which is proved by showing 



t sinh t 

< 



1 + 1 ~ 1 + sinh t 

for i > (put x = log(sinht)). Set p{t) = t/(l + 1). Since p(t) is monotonically increasing, we 
have the desired inequality: p(t) < p(sinht). Hence (|5.4p is proved. 

Next, consider (|5.3|) . which is proved by showing |g(C)| < V% and 1 1 /c/(C) I < V^j where 

arcsinh(e^) 1 + e^ 

9{C) ~ 



1 + arcsinh( e^) e^ 



The functions g and 1/g are analytic on & n / 2 (and continuous on S$ n /2)- Therefore, by the 
maximum modulus principle, \g{C)\ and |l/g(C)| have their maximum on the boundary of @ n /2, 
i.e., | Im£| = 7t/2. It is sufficient to consider z = x + i (jl/2) for 

Let us show l/|g(,z)| < a/2 first, which is relatively easy. Setting X = log [ 1 1 — e 2x \ + e 2x ] 
and Y = arg[\/l — e 2x + i e x ], we have 



| 5 (z)| 2 = | 5 (x + i7T/2)| 2 = (l+e- 2 ^ X +1 " 



(X + i)2 + y2- 



Furthermore, from 



flog[2e 2a; - 1] (ifx>0), 
\0 (if x < 0), 

J tt/2 (if x > 0), 

\arctan(e a; /v / l - e 2x ) (if x < 0), 



it holds for x > that 

1 1 (1 + log[2e 2:r - l]) 2 + (tt/2) 5 



| 5 (z)| 2 1 + e" 2 * log 2 [2e 2a; - 1] + (tt/2) 2 

< i (i + {vr+^-i}/2) 2 + (V2) 2 

-l +e -2x ({^IT^-l}/2) 2 + (7t/2) 2 



and it holds for x < that 



1 1 1 1 f 1 

'i + — , = — = < — — + 



\g(z)\ 2 1+ e" 2:c \ axctan^e^/Vl- e 2x ] j ~ 1+ e" 2 * I (e*) 2 

Thus |l/ff(x + i7t/2)| 2 < 2 holds for all x G R. 

Next, let us show |<?(z)| < \/2- It holds for x > that 

. , ,„ log 2 [2e 2:r - 11 + (tt/2) 2 

For x < 0, we have 

-(.! + (-■'■) 



2 _ 2x , arctan 2 [e a, /vi — e 2x ] 



1 + arctan 2 [ e x — e 2x ] ' 
and let us bound the right hand side in the following three cases: 

a) -ilog[l + ^] <x<0, 
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b) -\ log[l + < x < -\ log[§] (note that -\ log[l + < -± log[§]), 

c) x <-|log[l + ^]. 

In the case of a), it holds that 



| 5 (z)| 2 <(l+e- 2 *) {n/2) 



In the case of b) , it holds that 

\g(z)\ 2 <(l+e- 2 *) 



In the case of c) , it holds that 



1 + (tt/2) 2 
1 + (tt/3) 2 



2 < ri+ri+4 



{n/2f 



<(i + ii + H 



1 + (tt/2) 2 
(tt/3) 2 



2. 



7tV/ l + (7t/3) 2 



2. 



Uj(z)| 2 <(l+e- 2 *) [eX/Vl ' 



2^1 2 



1 + [ex/V^O] 2 



1 



tanhx 



-2( 1 tanh(-ilog[l + i|]) 



7t 

1 + — < 2. 
18 



Thus |<7(rr + i tt/2) | 2 < 2 for all This completes the proof. 



Lemma 5.5. For all £ G and i £ M, we have 



< 



V2 



e C + Vl+ e 2 ?| ~ |1+ eC| ' 



1 



< 



1 



e x + y j 1 + e 2x - l + ' 
Proof. Consider (|5.6|) first, which is proved by showing 

1+ e x 



g(x) 



e x + VT+~e 



3x 



< 1 



for x£l, Since 



g'(x) 



e x 1 



,2.c 



(1 + e x ) - VT + ^j < 



(5.5) 
(5.6) 



we have liiiL r _-)._ 00 {g(x)} = 1, which is the desired result. 

Next, consider (|5.5|) . which is proved by showing \g(()\ < v2- The function g is analytic on 
@ n /2, (and continuous on S> n /2)- Therefore, by the maximum modulus principle, \g(C)\ has its 
maximum on the boundary of & n /2, i.e., | Im (\ = n/2. It is sufficient to consider z = x + i (n/2) 
for x£l, In the case x > 0, we have 



i 



In the case x < 0, we have 



2d' 



i + Vi 



-2x 



1+ e" 



-2.1' 



(1 + vT 



-2x\2 



< 1 + e~ 2x < 2 



\g(z)\ 2 = y/1 - e 2x + e 2x - i (1 - VI 
This completes the proof. 



e 2x ) e x 



1 + e 2x < 2. 
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Lemma 5.6. For all £ G 5^/2, it holds that 

1 



< 



1 



|1+ eC| - |l+ e 2 qV2- 
Proof. Let rr, y £ M with |y| < 7t/2, and let £ = x + iy. Then we have 



1 




1 


1+ e 2 f 




1 + eC 



'1 + 2 cosh(a;) cos(y) 



2(cosh(x) + cos(y)) y cosh (x) — sin 2 (y) 



>0, 



which completes the proof. 

Lemma 5.7. Let x, y € M with |y| < 7t/2. Then it holds that 

cosh 2 [(2; + iy)/2] 



sup 



cosh(x + iy) 



< cos 2 (y/2) = 1 ^ | I 
cos y 2 



cos y 



Proof. Firstly, we have 

cosh 2 [(x +iy)/2] 



cosh(x + iy) 



cosh(x) + cos(y) 



< 



cosh(x) + cos(y) 



2 J cosh 2 (x) - sin 2 (y ) 2 J cosh 2 (x){l - sin 2 (y) } 



which then is estimated as 

cosh(x) + cos(y) 



cosh(x) + cos(y) 1 + cos(y) 



2^cosh 2 (x){l-sin 2 (y)} 2 cosh(x) cos(y) 2-l-cos(y)' 

This completes the proof. 

Then Lemma 15.31 is proved as follows. 
Proof. From the inequality (|2.6|) with i = 3, it follows that 



\m\<K 



arcsinh( e^ 



1 + arcsinh( e^) 



a-l 



1 



eC + Vl + e 2 ? 



1 



1 + e" 2 C 



1/2 



Consider the case 1 < a first. From Lemmas 15.41 15.51 and 15.61 it holds that 



1^(01 < K 



V2 



a-l 



1 + e~< 

< 2 (a+/3-l)/2 K 



1 + eC 



1 + e" 2 C 



1/2 



1 


(o-l)/2 


1 


0/2 


1 


1/2 


1 + e" 2 C 




1+ e 2 C 




1 + e" 2 C 





for all (Gfd. For x G R, it holds that 



\F(x)\ < K 
< K 



1 



1 + e~ x 
1 



a-l 



1 



1 



,1 + e~ 2x 
This completes the proof for 1 < a. 



1 + e x ) V 1 + e~ 2x 

i..-l)/2 / j x jfl/2 



1/2 



1 + e 2a; 



1+ e- 



-2x 



1/2 



15 



Next, consider the case < a < 1. From Lemmas 15.41 15.51 and 15.61 it holds that 



1^(01 < K 



K 



1 + e~< 



l-a 



V2 

[1 + e"9 2 



V2 



2(1 + e- 2 ? ) 



1 + eC 

(l-a)/2 



1 + e-2< 
^2 



1 + c< 



1/2 



1 



Furthermore, from Lemma 15.71 it holds for £ = x + iy 6 



(1+ e-C) 2 




cosh 2 (C/2) 


2(1 + e" 2 9 




cosh(C) 



1 



<^ 1 + 



1 + e" 2 ? 
that 

1 



a/2 



cosy 



Then using Lemma 15.61 we have 

\F{()\ < (l + 



cos d 



(l-a)/2 



1 



/3/2 



for all C G For x € R, it holds that 

l-a 



|F(x)| < K 



K 



1 + e -x 

' (1 + e- 
2(1 + e" 2a; ) 



1 



1 + e^ 

■x\2 \ (!- Q )/ 2 



1 



1+ e 2 ? 
1 

1 + e" 2 * 



1 + e" 2 < 



a/2 



<K[- 1 + 



cos 



(l-a)/2 



1 + e x 
1 



1/2 



1 



1 + e 2 * 



1 + e" 2:r 
0/2 / x 



a/2 



1+ e- 



-2x 



a/2 



This completes the proof for < a < 1. 



5.2 In the Case of the DE Transformation 

We need the following definition in the case of the DE transformation. 

Definition 5.2. Let L, R, a, f3 be positive constants, and d is a constant with < d < n/2. 
Then L£ E ^ a p(@d) denotes a family of functions F that are analytic on f^, and satisfy for all 
(Gfj and x£l that 



\F(0\ < 
\F(x)\ < 



(7T/2)L|coshC| 



|]_ _|_ g— nsinhf |a/2M _|_ grcsinhf |/3/2 ' 

(7t/2)i?cosha; 



(5.7) 
(5.8) 



^ _|_ g— 7tsinha;^a/2^ _|_ g7isinha;^/3/2 ' 

If F belongs to this function space, the errors of the Sine quadrature and the Sine indefinite 
integration are estimated as below. The proofs are omitted since they are quite similar to the 
existing theorems for the case 4 [3J Theorems 2.14 and 2.16]. 

Theorem 5.3. Let F G a p(®d) i let fj, = min{a, /?}, let v = max{a, /?}, let h be defined 
as (|2.9|) . and let M and N be defined as (|2.1U|) . Furthermore, let n be taken sufficiently large 
so that n > (ve)/(8d), Mh > x a , and Nh > xp hold. Then, 

AT 



/oo 
F(x)dx-h F ( kh ) 
-°° k=-M 



< 



/i 



2L 



(1 - e-^ e / 4 ){cos(f cos d)}^)/ 2 cos d 



+ #e™/ 4 



-2Ttdn/ log(8dn//i) 
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Theorem 5.4. Let F £ L^ija pi^d), let /x = min{a, /?}, let 1/ = max{a, let /i be defined 
as (|2.11|) . and let M and N be defined as (|2.1U|) . Furthermore, let n be taken sufficiently large 
so that n > (Ve)/(4cf), M/i > x a , and iV/i > xp hold. Then, 



sup 



^ N 

/ F(x)dx- 51 F{kh)J(k,h){£) 



k=-M 



< 



fid 



L 



(1 - e-^ e / 2 ){cos(f sin d)}^)/ 2 cos d 



+ Re n(a+fS)/4 



In view of Theorems 15.31 and 15.41 our project is done by checking F G L£ E ^ a s(@d) i n each 
case: 1, 2, 3. Let us check them one by one in what follows. The next lemma is useful for the 
proofs. 

Lemma 5.8 (Okayama et al. Lemma 4.22]). Let x, y € R with \y\ < n/2, and let C = x + iy. 
Then, 



1 



^ _|_ g7tsinh f 
1 



^ _|_ g— nsinh^ 



< 



< 



(1+ e 7tsinh ( :r ) cos y)cos(f siny)' 
1 

(1+ e- nsinh ( x ) cosy) cos(% sin y)' 



5.2.1 Proofs in the Case 1 (Theorems [277] and [33]) 

The claims of Theorems 12.71 and 13.31 follow from the next lemma. 

Lemma 5.9. Let the assumptions in Theorem 12 .71 or Theorem l3.3l be fulfilled. Then the function 
F(C) = /OWCMeiCC) belongs to L£ e Rq/3 (^) with L = 2^/{cos(f sin d)}^V 2 and 
R = 2 V K. 



Proof. Consider the case a < f3 first. From the inequality (|2.ip . it follows that 

K 



\F{Q\ < 



2 a l + e 



|]__|_ g— 7tsinhf |q/2 | ]^ _|_ g7tsinhf |/3/2 

for C € ^ with Re£ < 0, and from the inequality (|2.2|) . it follows that 

K 2? 

\F{Q\ < 



7tsinhC|(/3-a)/2 



|]__|_ g— 7tsinhf |o/2 1 ^ _|_ g7Tsinhf |/3/2 1 1 + g — ^ s i n h C | a )/2 

for ( 6 ^ with Re£ > 0. Setting £ = x + iy with x < 0, we have 
2 a \l 4- g^sinhfK/ 3 -")/ 2 



2 Q (1 + e 



7Tsinh(x) cos i/\ (f3—a)/2 



sin 2 (7t cosh (x) siny) 
cosh 2 (7rsinh(x) cos y) 



< 2 Q (1 + e °) (/3 - a)/2 {1 - 0} (/3 ~ q)/4 = 2 (a+/3)/2 < 2^ < 



08-a)/4 



2^ 



{cos^siny)}^"")/ 2. 
Furthermore, setting £ = x + iy with x > and using Lemma l5.8t we have 
2^ 2^ 



< 



|1+ e -7tsinhC|(/3-a)/2 _ (1 + e -7rsinh(x) cos y ) (/3-a)/2 { COS (| sin y)}^-")/ 2 

2? 



< 



(l+0){cos(f siny)}^"")/ 2 
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Thus, since \x = a and v = j3 in this case, it holds for all C, £ S> d that 

if 2^ 



\F{Q\ < 



1 1 + e-« sinh C|«/2|i + e " sinh C | P/2 {cos ( f sin d) } ("-/0/a ' 
and it holds for all x S M that 

K 2 V 



\F(x)\ < 



(1 + + e nsmhx^/2 { COS (f sin 0)}>^)/ 2 ' 

In the case a > f3, the same inequalities hold. This completes the proof. ■ 

5.2.2 Proofs in the Case 2 (Theorems [278] and ETift 

The claims of Theorems 12.81 and 13.41 follow from the next lemma. 

Lemma 5.10. Let the assumptions in Theorem 12.81 or Theorem 13.41 be fulfilled. Then the 
function F(() = /(^ DE2 (CM E2 (C) belongs to Lg^ ^) with L = K and R = K. 

Proof. From the inequality (|2.6p with i = 2, immediately (|5.7|) and (|5.8p hold with L = R = 
K. U 

5.2.3 Proofs in the Case 3 (Theorems \2J)\ and [33]) 

The claims of Theorems 12.91 and 13.51 follow from the next lemma. 



Lemma 5.11. Let the assumptions in Theorem 12.91 or Theorem 13.51 be fulfilled. Then the 
function F(C) = /WWOVWO belongs to l^\ 2a ^{9 d ) with L = 2{c d ) l - a K and R = 
2(e n / 12 ) 1 ' a K, where c d is a constant defined in (|2.12|) . 

For the proof, we need the next estimate. 

Lemma 5.12. Let d be a constant with < d < n/2. Then, 

1 + log(l + e nsinh<: ) 1 



sup 



sup 



log(l + e nsinh ^) 1 + e~ nsinh ( 
1 + log(l + e 718 ^) 1 ] 



logfl + e 71 s ™' 1 x ) 1 + e — 71 s ™' 1 x 



< c d , (5.9) 

< e n/12 , (5.10) 



where c<i is a constant defined in (|2,12[) , 

Proof. Consider (|5.1Up first, which is proved by showing 

1 + t n -t\^i lo S 6 



! ,( t ) = i±I(l-e-')<l + ^ (<e»/-) 



t v ' - 6 

for t > (put t = log(l + e nsinhx )). Let A be a value with log 6 < A < log 7 so that p'(X) = 
(1 + A + A 2 — e A )/(e A A 2 ) = 0. Then clearly p(t) < p(X) holds. Furthermore, using the relation 
1 + A + A 2 = e A , we have 



Since the function 1 + ie x is monotonically decreasing, 1 + Ae A < 1 + (log6)e log6 holds 
(note that log 6 < A). Hence (|5.10p is proved. 
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Next, consider (|5.9|) . By the maximum modulus principle, that is proved by showing 



where £ = log(l + e 7rsmh ( a; + lc 0) an d j = — log [cos sin d)]. Here, notice that 

Re£ = log |1 + e nsi,oh ^ +i ^ \ > log [(1 + e ™inh(a)cosd) cos (f sind)] > -7 

holds from Lemma [5,81 Let us bound |p(£)| in the two cases: a) |£| < log(2 + e 7 ), and b) 
|£| > log (2 + e 7 ). In the case of a), it holds that 



\p(0\ 



(1+0E 



(-0 



fc-i 



k=l 



kl 



<(i+i^i)f: K ^ = 1 |^(e |?i -i)- 



k=l 



Furthermore, since (1 + x)( e x — l)/x is monotonically increasing, we have (|5.1ip . In the case of 
b), since Re£ > 7, it holds that 

\P(0\ < ^(1 + I e-«|) = ^(1 + e-^) < i + + e 7 ). 

Furthermore, since (1 + x)/x is monotonically decreasing, we have (|5,11|) , This completes the 
proof. ■ 

Then Lemma 15.111 is proved as follows. 

Proof. From the inequality (|2.6|) with i = 3, it follows that 



1^(01 < 



1 + log(l + e nsinh<: ) 



1-a 



log(l + e nsinh <) 
l+log(l+ e nsinh<: 



|1 -)- gTtsinhf |/3 
1 



log(l + e ns ' mh () 1 + e - ns ' mh ( 
Then use Lemma 15. 121 to obtain the desired result. 



K(7r/2)|coshC| 



I ^ I g— Ttsinhf laM _j_ g7tsinhO/3° 
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